Device-independent bounds for Hardy's experiment 
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In this Letter we compute an analogue of Tsirelson's bound for Hardy's test of nonlocality, that 
is, the maximum violation of locality constraints allowed by the quantum formalism, irrespective of 
the dimension of the system. The value is found to be the same as the one achievable already with 
two-qubit systems, and we show that only a very specific class of states can lead to such maximal 
value, thus highlighting Hardy's test as a device-independent self-test protocol for such states. By 
considering realistic constraints in Hardy's test, we also compute device-independent upper bounds 
on this violation and show that these bounds are saturated by two-qubit systems, thus showing that 
there is no advantage in using higher-dimensional systems in experimental implementations of such 
test. 



Introduction.- The development of quantum infor- 
mation science is based on a recurrent pattern: non- 
classical features of quantum physics, previously con- 
sidered as mind-boggling and worth only of philosoph- 
ical chat, are found to have an operational meaning 
and even to be potentially useful for applications. One 
of the discoveries that triggered this development is 
the prediction and observation of the violation of Bell 
inequalities This observation implies that correla- 
tions obtained by measuring separated quantum sys- 
tems locally cannot be simulated classically without 
communication, a fact that is often referred to as non- 
locality. 

Within quantum information, nonlocality has under- 
gone an interesting parable. For many years, it has 
been put aside as having fulfilled its role: the loathed 
local variables models having been disposed of forever, 
one could peacefully concentrate on entanglement the- 
ory. Only few researchers kept on believing that this 
very intriguing observation could be useful for some- 
thing in itself. The latter view was vindicated a few 
years ago, when it was noticed that nonlocality allows 
device-independent assessments: indeed, nonlocality is 
assessed only from the input-output statistics of the 
measurement, without reference to the degree of free- 
dom that is being measured. This powerful type of as- 
sessment is sensitive to the existence of undesired side- 
channels and will be ideal for certification of future 
quantum devices. So far, device-independent results 
are available for the security ofguantum cryptography 
[H, 13] , the quality of sources [J Q and measurement 
devices (gL the amount of randomness that one can 
generate [3, Q . In this paper, we study the possibility 
of device-independent assessment of one of the earli- 
est proposals to check nonlocality: it used to be called 
Hardy's paradox but, in the spirit of quantum informa- 



tion, we'd rather call it Hardy's test [g[ 

Hardy's test was originally stated by means of a par- 
ticular experimental setup consisting of two overlap- 
ping Mach-Zehnder interferometers, one for electrons 
and one for positrons, arranged so that if the positron 
and the electron each take a particular path they will 
meet and annihilate one another. A paradox arises 
under the assumption of local realism: in any classi- 
cal local theory a certain detection pattern must never 
occur, while quantum theory assign to its occurrence 
a nonzero probability, hereafter referred to as Hardy's 
probability. It was soon realized that the argument 
could be extended to different states and measurements 
(ill 11 1, and proved to hold for almost all entangled 
pure states of two qubits, with maximum Hardy's prob- 
ability equal to 



PHardy 



= (5VE-II) /2 ps 9% 



(1) 



Interestingly, though, the maximaly entangled state of 
two qubits does not show nonlocality in Hardy's test. 

Hardy's test has been the object of several theoret- 
ical generalizations 12-2^] and has been im plem ented 
in experiments using photonic systems 2ll425|. The 
latter, however, had to consider deviations from the 
original proposal, where the probabilities of a set of 
observations - hereafter referred to as constraint prob- 
abilities - were assumed to be strictly equal to zero, an 
obviously unrealistic requirement. One way to over- 
come this problem is to consider a nonidcal version of 
Hardy's test, and to compute local bounds on Hardy's 
probability in terms of relaxed bounds on the con- 
straint probabilities. The computed local bound, which 



a successful experiment must violate, turns out j26H28| 
to be equivalent to the Clauser-Horne (CH) Bell in- 
equality [III. 

In this paper, we provide three device-independent 
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FIG. 1: Schematic diagram for the Hardy's test scenario. 

results on Hardy's test. First, we consider the original 
(or ideal) Hardy's test and prove that that Q) is the 
maximum value of pHardy allowed by quantum physics, 
irrespective of the dimension; this is the analog of the 
Tsirelson bound [3(| ■ A remarkable consequence of 
our derivation constitutes our second main result: any 
state that achieves ([T]) in the ideal test is equivalent, 
up to local isometries, to the unique two-qubit state 
that achieves that violation. This is a case of self- 
testing 0, [3l| , the first that detects a non- maximally 
entangled state (see parallel work [HJ]). Finally, our 
third result is a proof that, even for nonideal versions 
of Hardy's test, there is no practical advantage in using 
higher-dimensional systems. 

Hardy's test -Let us briefly summarize Hardy's test. 
Consider two parties, say, Alice and Bob, each of which 
is able to perform two possible measurements, x — 
{Aq, A\} and y = {Bq, B\}, respectively, on its part of 
a shared physical system. Each measurement has two 
mutually exclusive outcomes, labeled by a = {±1}, for 
the measurements of Alice, and b = {±1}, for the ones 
of Bob. The situation considered by Hardy assumes 
the three constraint probabilities 



P {+,+\Aq,B ) = 0, 
p{+,-\A 1 ,B Q ) = 0, 
p(-,+\A o ,B 1 ) = 0. 



(2a) 
(2b) 
(2c) 



Suppose there are measurement devices and physical 
systems such that these three equations are fulfilled. If 
this setup can be described by a local realistic theory, 
then it follows that 



PHardy = p(+,+\ A l> B l) = °- 



(3) 



Hardy realized that in quantum mechanics there are 
measurements and a particular state of a two-qubit sys- 
tem such that the constraint probabilities are fulfilled 
while Hardy's probability is nonzero, leading to a so- 
called 'paradox'. Extending the analisys to all possible 
measurements and states, Hardy later showed that the 
maximum value of pHardy for systems of two qubits 
is (5\/5 — ll) /2. A brute force calculation proved 
that this value cannot be exceeded using two three- 
dimensional systems Here, we prove that this 



value is device-independent, that is, it is optimal for 
bipartite quantum systems of any dimension. 

Theorem 1. The maximum value of Hardy's probabil- 
ity for quantum systems of arbitrary finite dimension 
is PHardy = (5\/5 — ll) /2, just as for qubits. 

Proof. In quantum mechanics, joint probabilities for 
the outcomes of measurements performed on space- like 
separated parts of a quantum system are given by 



p(a, b\x, y) = Tr (pU a \ x ® n 



b\v) 



(4) 



where p is the state of the system and Ha^n^ are 
the measurement operators associated to outcomes a, b 
of measurements x, y, respectively. The latter opera- 
tors are POVM effects, in general; however, since we do 
not set any constraint on the dimension of the Hilbert 
space, Neumark's theorem allows us to consider only 
projective measurements, without loss of generality. 
The core of the proof exploits the following lemma, 
proven in 



33 



Lemma 1. Given two Hermitian operators Aq and A\ 
with eigenvalues ±1 acting on a Hilbert space %, there 
is a decomposition of % as a direct sum of subspaces 
W" of dimension d < 2 each, such that both Aq and 
A\ act within each H 1 , that is, they can be written as 
Aq = ® t A and Ax = (& i A[, where A and A\ act 
onW. 



Let then Aq = H_ 



\A 



U-\A and A x = H + \ Al 



n_| J 4 1 , where H a \ x arc projection operators. It follows 
from Lemma 1 that = 0^^, where each 
acts on TV , for all a and x; we also denote IF = Ii' l +1 ^ x + 
rP-i|a: the projector on H l . Needless to say, Lemma 1 
is also valid on Bob's side; we use analog notations for 
Bob's operators. With these notations, 

p(a, b\x, V)=J2 ^ Tr (phK\x ® U l\ y ) ( 5a ) 

= ^2qtjPtj( a i b \ x ,y)i ( 5b ) 

where qij ~ Tr (pTF CED IP') and pij = 
(IF (g> H'jon 1 (g) IP') /qij is, at most, a two-qubit 
state. Since q^ > for all i,j and YlijQij = 1; the 
constraint probabilities $2$ are satisfied for p if and 
only if they are satisfied for each of the py. But, then, 

p(+,+\A 1 ,B 1 )=J2<lijPij(+,+\ A i, B i), (6) 

is a convex sum of Hardy's probabilities in each two- 
qubit subspace [Hf . As a convex sum, it is less or 
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equal to the largest element in the combination, whose 
maximum value is known to be given by ([1}. This 
concludes the proof [39| . □ 

Hardy's test leads to self-testing. - It follows from 
the previous proof that p(+, +\Ai, B{) reaches its max- 
imal value if and only if pij(+, +\Ai, B\) is maximal 
for every ij such that qij 7^ 0. The following Lemma, 
proved in [13, [Hi , states that only a very specific class 
of two-qubit states can lead to this maximal value: 

Lemma 2. Consider Hardy's test implemented in a 
two-qubit system, and let Aq = Bq = |0) (0| — |1) (1|. 
The probability pHardy reaches its maximal value if, and 
only if, the state of the system is 



a (|01) + 1 10)) + e'Vl - 2a 2 |11) 



(7) 



the 
Bi 



other 



two 



and 
A t 

^(VT=2P|0)-e"a|l)), 
and 9 is arbitrary. 



measurements are 
-){-| wi th |+) = 

(3 - V5) /2 



In view of this, one can conjecture that, if the max- 
imal value of PHardy is observed, the state must some- 
how be a direct sum of copies of \<j>). We proceed to 
prove that this is indeed the case: 

Theorem 2. If p Hardy = (5v5 — ll)/2 is observed in 
an ideal Hardy's test [i.e., together with ©/, then the 
state of the system is equivalent up to local isometries 
t° \°) ab ® l^A'B'' w here \<j>) is given in Qanrf \a) is 
an arbitrary bipartite state. In other words, the ideal 
Hardy's test constitutes a self-testing of \<f>) . 

Proof. Without loss of generality, let us choose the 
eigcnbases of Aq and Bq as the computational bases: 



n 



+ |A 



120 (2il n 



\2i + i) 2i + i|, n J 



2j + 1) (2j + Then, by Lemma 



\2j) (2j\, ILl ]Bo 
m JHj(+,+\^l,Bi) = Tr ( PlJ W MAi ® n J +|Bi 



(5\/5 — ll) /2 if and only if pij = \<j>%j) (< 
\ ( f> ij ) = a(\2i,2j + l) + \2i + l,2j)) + 



where 



e i Vl-2a 2 |2i + l,2j + l), (8) 



and a = y (3 — y/S) [2 and arbitrary 0. This way, a 

state \tp) can lead to a maximal value of pHardy if, and 
only if, it is given by 



(9) 



The coefficients qij are arbitrary probabilities that, by 
definition, are constrained to the form q^ = riSj, where 
r i,Sj > 0, r « — s j — 1- The angle cannot 
depend on the indices i,j because IFm^ is uniquely 
defined by 9 (cf. Lemma [5]), and, by definition, is 
independent of j; the same reasoning can be applied 
to Now, following [f|, we append local ancilla 

qubits prepared in the state |00) A , S , and look for local 
isometries $a and $ b such that 



{* A ® * fl ) |V) AB |00) 



I A'B' 



I AB 



(10) 



where |cr) is a bipartite 'junk' state. This can indeed 
be achieved for <&a = $5 = $ defined by the map 

<f|2fc,0) cc , »\2k,0) cc ,, (11a) 
$|2fc + l,0) cc , ^|2fc,l) cc ,, (lib) 

for both C = A,B. □ 

Up to now, self-testing was known only for max- 
imally entangled states (see, e.g., [f| and references 
therein). A parallel, independent work by Yang and 
Navascues provides a very general approach to the 
self-testing of bipartite non-maximally entangled states 
(32I ] . Remarkably, though, our Hardy point is not de- 
tected by that test j40T |. 

Hardy's experiment with realistic constraints. - Sup- 
pose now that the constraint probabilities ([2]) in 
Hardy's experiment are not exactly equal to zero. In 
this case, the local bound on Hardy's probability is no 
longer zero, either: in general, it is given by the follow- 
ing inequality 



p(+, +\A lt Si) < p(+, +\A , B )+ 

p(+,-\A 1 ,B Q )+p(-,+\A ,B 1 ). (12) 

This inequality is a re-writing of the CH inequality 
[29j], which is not surprising, since the CH inequality 
is the only relevant criterion for nonlocality in a sce- 
nario with two parties, two inputs and two outcomes. 
In other words, as noticed in |28j . Hardy's experiment 
turns out to be a study of the violation of the CH in- 
equality under further constraints about the values of 
some probabilities. 
Let us now set 



p(+,+\A ,B ) < e, 
p(+,-|i4i,B ) < e, 
p{-,+\A ,B 1 )<e, 



(13a) 
(13b) 
(13c) 



for some e > [41|. The local bound on Hardy's prob- 
ability becomes 



L4i,Bi)<3e. 



(14) 
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For e > |, the bound is trivial and quantum physics 
certainly cannot violate it; while for < e < |, quan- 
tum physics may lead to a violation of the local bound. 
As before, we want to assess the maximal quantum 
violation in a device- independent scenario, i.e., with- 
out making any assumption on the Hilbert space di- 
mension. The previously stated theorem cannot be 
extended, so we take a different approach: first, we 
use semi-definite programs to obtain an upper bound 
on Hardy's probability, using the method of Navascues, 
Pironio and Acin [34| ; second, by considering two-qubit 
systems we obtain a value that is certainly achievable 
with quantum systems. By noticing that the values 
thus obtained coincide, we conclude that we have ob- 
tained the optimal value for Hardy's probability, and 
that this value can be reached with two qubit systems. 

In detail: let Q be the set of quantum joint proba- 
bility distributions, that is, vectors of probabilities of 
the form (TJJ, for all a,b,x,y. We compute a device- 
independent upper bound on Hardy's probability by 
optimizing it not over quantum probabilities in the set 
Q but over a larger set of probabilities that is compu- 
tationally tractable — as opposed to Q, that still lacks 
a better characterization. This set is one of an infinite 
hierarchy of sets Qi D Q2 D ■ ■ ■ D On D ■ ■ ■ , defined 
in terms of semi-definite programs 34, 35|, proven to 
converge to the quantum set, lim„_>. 00 Q„ = Q. For 
several values of e in the interval 0<e<l/3, we op- 
timize Hardy's probability over the set Q3, enforcing 
the constraints (|13[) . The implementation was done 
in MATLAB using semi-definite programming 
The results form the solid line in Fig. [2j For the lower 
bound, we consider the most general mixed states of 
two qubits and POVM elements acting on those. The 
maximal value of the Hardy's probability is estimated 
using constrained nonlinear optimization methods in 
MATLAB. These methods are not guaranteed to con- 
verge to global maxima, though, and are in fact rather 
sensitive to seed conditions; each point on the dotted 
line in Fig. [2] is the maximum obtained over 10 4 runs, 
with random initial seeds. 

The computed lower and upper bounds for Hardy's 
probability differ, at most, by values of order 10~ 2 ; in 
the region e < 0.2 (where any experiment that aims 
at implementing Hardy's test will have to be), this 
difference is of order 10~ 6 . This proves that there is 
no advantage in using higher-dimensional systems, as 
compared to two-qubit systems, even in the presence 
of imperfections. 

Conclusion.- In this letter, we prove that the maxi- 
mum value of Hardy's probability found for two-qubit 
systems, (5\/5 — ll) /2, is the maximum one allowed 



FIG. 2: Upper and lower bounds on maximum Hardy's 
probability pHardy in terms of the bound e on the con- 
straint probabilities. The solid (blue) line is the upper 
bound, computed from the set Q3; the dotted (red) line 
is the lower bound, computed from two-qubit systems; the 
dashed (black) line is the local bound. 



by quantum theory, irrespective of the dimension of the 
system and of the measurements performed, that is, in- 
dependend of the devices used. By showing that only a 
certain class of states can lead to such maximal value, 
we show that Hardy's test is, in fact, a self-testing pro- 
tocol for such states. Extending the first results to a 
nonidcal vesion of Hardy's test, where the constraint 
probabilities are no longer equal to zero, we compute 
device- independent upper bounds on Hardy's probabil- 
ity, in terms of the error parameter, and show that this 
bound is saturated by two-qubit systems. 

Despite their fundamental importance, as the first 
proven analogue of Tsirclson's bound for Hardy's test, 
the results here presented also serve as a guideline 
for future experimental implementations, as they show 
that there is no advantage in using higher dimensional 
systems, as compared to two-qubit systems. 
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It is important to note that for the maximum value of 
([6]) to be reached it is necessary that, for all i,j such 
that qij 7^ 0, the dimension of both T-C and W be equal 
to 2. This implies that the effective dimension d of the 
local Hilbert space H of the system is even, and that 
Ao, Ai, Bo, B\ have exactly d/2 positive and negative 
eigenvalues. 

There is an alternative, simpler proof of the above the- 
orem that consists, basically, in noticing that any prob- 
ability distribution that maximizes Hardy's probabil- 
ity is an extremal point of the set of quantum proba- 
bility distributions. According to [3^], every extremal 
point, in this scenario, can be obtained from projective 
measurements on two-qubit systems, thus proving the 
stated result. The reason we opted for presenting the 
extensive proof is that it leads to interesting insights 
about the states that lead to such maximal violation, 
as discussed below. This proof cannot be extended to 
the nonideal scenario we later consider due to the fact 
that, in that scenario, it is not clear wether the points 
reaching maximal Hardy's probability are extremal or 
not. 

At least not up to the Q2 step of the hierarchy (T.H. 
Yang, private communication, 17 July 2012) 
Notice that, if no-signaling holds, then 
p(+,+\A ,B ) = e implies p(+\A ) > e and 
p( — , +|^4o, B\) = e implies p(— \Ao) > e. Therefore 
e < 2- But the region of interest is in fact e < | as 
explained just below. 



